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£h ' Abstract 
3. 

Keeping the discreteness of the reciprocal space we calculate the spectrum of incoherent electron-hole excitations 
in the conducting Fermi liquids. The method is illustrated on the well-known jellium model within the random phase 
approximation. It also leads to the formulation of a sum rule from which we get the details of the dispersion curve 
for the collective plasmon mode. The notion of time averaging in the discrete approach is briefly recalled. 
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(ZJ ! I. INTRODUCTION 

The electron-hole excitations in conductors invoke in the frequency dependence of dielectric function e(k, u>) a 
dense alternation of poles and zeros at the scale of discreteness of the reciprocal space Usually this complex 
■"^J , mathematical problem is avoided by making the standard procedure of continuation of the wave vector variable k 
£h ' which is well founded for macroscopic systems, and in addition includes the standard proposition allowing for the 
O \ irreversibility in the thermodynamic limit ||. By this the dense set of poles and zeros is eliminated and e(k, uj) 
becomes an analytic functions. In particular, within random phase approximation (RPA) for jellium which will be 
followed here, Im e is a continuous function of k and uj in the range of variables covering the so-called electron-hole 
continuum. 

Although cumbersome at first sight, the discrete presentation of dielectric response still appears to be advantageous 
\q ' and physically more transparent in some treatments, like in the calculations of cohesive energy j^^J and one-particle 
£NJ , spectral function lq] in the single band and multiband systems. In the present work we show how the method 
^s0 developed in Refs. reproduces some well-known results for the jellium model. In particular we derive an explicit 
^ \ expression for energies of incoherent excitations, which is of course not attainable after continuation. This leads to 
the formulation of a sum rule which enables a direct determination of the dispersion curve for the collective excitation 
(i. e. plasmon in the present example). The particular detail which then can be followed in a transparent way is the 
d ' cross-over of the plasmon dispersion from the collective to the incoherent regime as the wave vector increases, the 
£h ! subject which was often exposed incorrectly or imprecisely in literature. 

1 ■ In Section II we start with the formulation of the problem and continue with the explicit calculation of the energies 

of incoherent excitations. Section III contains the formulation of the sum rule, the analysis of the plasmon dispersion 
g ■ and the short note on the calculation of Im e in the discrete approach. Concluding remarks are given in Section IV. 

O ■ 

II. INCOHERENT EXCITATIONS 

X: 

We start from the well-known RPA dielectric function for the jellium |1J. The excitation energy for a given value of 
the wave vector q is the solution of the equation e(q, uj) = in the frequency variable uj. Let us write this equation 
in the form 

e(q,uj) = l — > g(ku,q)- r — — , , ' -^rr. r = 0, (1) 
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where L 3 is the volume of the system. We keep throughout this text the discrete summation in terms of the wave 
vector component parallel to q, 
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The equation ku = constant defines the locus of the constant electron-hole excitation energy 
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i.e. the summation in Eq.(Q) goes over all such loci. As shown in Fig.|T] the locus for given values of q and ku is 
either circular [for E(k\\,q) > (kp — fcjj)/2m] or annular [for E(k\\,q) < (kp — k 2 )/2m] surface, centered at k\\. The 

k?,, while the latter is bounded by concentric circles with radii k±M and 
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The number of the points multiplied with the elementary surface (2n/L) of the reciprocal space in a given locus is 
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The dielectric function (Q) diverges at energies which are infinitesimally close to the values of electron-hole excitation 
energies (^) at the real w-axis. The energy difference between neighboring poles (|3|) is 



AE(q) = E{k\\ + 2tt/L, q) - E(k\\ , q) 
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Between each such pair of neighboring poles (S) there should be a zero of dielectric function. In other words all 
solutions of the equation e(q, uS) = except the largest one lie between neighboring electron-hole energies. Let us 
denote these zeros by f2(fcii,q), and write 



n(fc|| , q) + isignCO)^ = E(k {l , q) + 6(fc|| , q)AE{q) 



(7) 



with < 0(fc||,q) < 1. Let us remind that in contrast to these incoherent excitations, the excitation with largest 
energy can be macroscopically (with respect to 1/L) far from its closest electron-hole counterpart, i?(/c_F,q). Such 
isolated zero is a collective excitation, i. e. plasmon in the present jellium model. 

An approximate expression for the energies of incoherent excitations follows after retaining the ^-dependence only 
in those terms of Eq.(Q) which contain nearest neighboring poles to a given zero 0(fcii,q), i. e. only in terms with 
E{k\\ , q) and E(kn + 2ir/L, q). In terms with wave vectors different from kn or fcii + 2ir/L we substitute ui + isign(ui)n 
with E(kii,q). The validity of this approximative step is discussed in the Appendix. Eq.(Q) now reduces to 
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where ciq — I /me 2 is Bohr radius and 
F(fc||,q 
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With the most divergent terms singled out, we can now make continuation of the residual sums (|9|), by replacing 
each term with the integral on the interval 2tt/L. Choosing to replace a term characterized by the discrete value fc|, 
with the integral from fcf, — 2ir/L to fc[, we come to the expression 
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where ku — 2it/L and fcii + 2tt / ' L are the lower and upper boundaries of the interval that corresponds to the missing 
terms k\\ and fe|| +2ir/L in the sum (|^). We note that the final result of integration in Eq. (|To|) does not depend on 
the choice of integration boundaries for 2it /L intervals, as specified above. 
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Since the primitive function in the above integral is odd with respect to the pole at , the expression ( fL0[ ) can be 
rewritten in the form 



g(kLq)E(kLq) 
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After elementary integration we get 



L.2 _ 7,2 
/i p /in 

-2fci? - — Mn 



(k F - k\\)(k F + fc|| + q) 



(k F + h\)(k f 



q) 



(2fcn + q) In 



(12) 



Furthermore, by solving the biquadratic equation (||) we get for the function 0(fc||,q) in Eq. (]?]) the expression 

@(^|| , q) = i[l + s-sign(a;)v / TT^], 



where 
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Evidently O varies in the range (0, 1) in accordance with the starting assumption. More precisely 

< e(fe| h q) < 1/2 >0 
9(fc|| , q) = 1 for n 2 q 2 a - 4F(fc|| , q) = . 

1/2 < 6(fc||,q) < 1 <0 
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The above results can be written in a more explicit way in the limits q << k F and q » 2k F . In the former case 
and for fen <C kp the expression (113) reduces to 
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and the leading term in the energy of incoherent excitations is 
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i.e. the zeros of e(u>, k) in the w-plane lie very close to the corresponding electron-hole excitation poles. Here and 
further on we omit for simplicity the infinitesimal imaginary part in excitation energies (0). 

As the wave number k\\ increases the energy of incoherent excitations moves gradually towards the first higher 
neighboring electron-hole excitation. For k\\ close to k F and for kp 3> q ^> k F — ku the function (n2\j reduces to 



n(2k F + q) 2kp + q nkp irk F (k F — k\\ 
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Keeping only the leading term in the above expression the approximate expression for corresponding energies of 
incoherent excitation reads 



n(k h q) = E(k\\ + 27r/i,q) 
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mL \a(2k F / q) 



(19) 



i. e. these energies are very close to the energies of the first larger neighboring electron-hole excitations. 
In the latter case q » 2k F Eq.(|i~2"|) reduces to 
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Energies of incoherent excitations are then given by 



3 



E(k {] ,d), (21) 

i.e. they are very close to their electron-hole counterparts for all values of fcii. The physical reason for this small 
departure of incoherent excitations from bare electron-hole ones is the weakness of the bare Coulomb interaction 
Aire 2 jq 2 in this regime. 

It follows from the foregoing analysis that the regime of small values of q qualitatively differs from that of large 
ones. In the former case the energy of incoherent excitations moves from their electron-hole counterpart towards the 
first larger neighboring electron- hole excitation, crossing the half- width A£'(q)/2 (Eq. as fc|| moves from the lower 
bound, — q, to the upper bound, kp — 2ir/L. On contrary, in the latter case it does not cross the half-width. It is 
interesting to estimate the value of the wave vector q cr which roughly separates these two regimes. To this end we note 
that in the regime of small values of q the difference between the energy of incoherent excitation and its electron- hole 
counterpart has the largest value for k\\ k F . Let us define the critical wave number q cr for which Q(k F ,q cr ) jl3| ) is 
equal to 1/2. In this case x in Eq. ( |l4j ) diverges, i.e. q cr is the solution of the equation 

Trq 2 cr a = -2k F + (2k F + q cr ) In 2kp + Qcr . (22) 

Qcr 

To summarize, for q > q cr the energies of incoherent excitations are closer to their electron-hole counterparts than 
to the next neighboring electron-hole excitations, regardless to the value of the wave number k». For q < q cr the 
energies of incoherent excitations are closer either to their electron-hole counterparts or to the first larger neighboring 
electron-hole excitations for small and large values of the wave number fen respectively. The wave number q cr roughly 
separates these two regimes. 
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III. COLLECTIVE EXCITATIONS 



In the above approximate calculations of energies of incoherent excitations we explore the fact that corresponding 
zeros and poles of the dielectric function densely alternate on the w-axis. This method however cannot be used for the 
excitation with the highest energy, lying above the maximum of electron- hole excitation energy E(k Fl q) for given q. 

In order to determine the Q(k F , q)-dependence of this collective mode, we write Eq. (|l|) in the form 

Et;=- kF 9Mm><i)ii K * kll l^-E 2 ^,^] -{q 2 ^/^ 2 )^ [u J 2 -E 2 (k h c l )] 

= 0. (23) 



Factors that multiply the highest powers, oj k F Lw and u> p 2 in the polynomial presentation of the nominator on 
the left-hand side are 

a kF = -q 2 Ln/4e 2 , (24) 

and 

a kF -2,/L= 9{k\\,<l)E(kh<l) + ^ L J2 ^falLl) ( 25 ) 

/t'l|— — kp — — kp 

respectively. As it follows from the elementary algebra the ratio ak F -2-n I hi a k F is equal to the sum of all zeroes of the 
Eq. (p3|). From this we get 

kp-2-n/L 

n 2 (k ll ,< i )-E 2 (k ll ,(i) + n 2 (k F ,c l )~E 2 (k F ,<i)=uj 2 pl . (26) 

k=—k F 



Above relation is a sum rule which states that the sum of differences of the squares of excitation energies and 
corresponding electron-hole energies equals to the square of the plasmon energy w 2 ; = AirNe 2 /{LPm) Q. The disper- 
sion of the collective mode £l(k F ,q) = fi p j(q) follows directly from this relation once the energies of the incoherent 
excitations are determined, as it was done in the previous Section. To illustrate, let us recall two characteristic points 
of the dispersion curve f2 p ;(q). First, since all incoherent and electron-hole excitations vanish for q = 0, the above 
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sum rule reproduces the well-known result that w p i is the energy of the highest long-wavelength excitation, i. e. of the 
plasmon Q. Second, it is known from the continuum approach |(| that the dispersion curve fl p i(q) touches the border 
of the electron-hole continuum at a finite wave number. Here we point out that this wave number coincides with 
the wave number q cr given by Eq. (p2f). The result by which plasmons do not exist as collective excitations at wave 
numbers at which all incoherent excitations are closer to their electron-hole counterparts than to next neighboring 
electron- hole excitations can be derived in the following way. The plasmon dispersion Sl p /(q) follows also from Eq. (||) 
after omitting the term with poles and substituting E(k\\, q) by f2 p ;(q) in the expression ([)]) for the function F(k\\, q). 
Thus we get 
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n 2 q 2 a 
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(27) 



In the continuum approximation the plasmon dispersion touches the border of the electron-hole quasi-continuum for 
^pi(ci) = (2k F + q)q/2m. By this Eq. ( f27j) reduces to Eq. j2^), i. e. the plasmon indeed ceases to exist as a collective 
excitation just at q = q cr . 

In the original works ^ as well as in the textbooks Jl|,||,0-||] the touching point q cr is usually interpreted as the 
wave number above which the decay of plasmons into electron- hole pairs takes part. Sometimes it is even stated or 
hinted that the dispersion curve fi p i(q) and the upper border of electron-hole range cross at q = q cr We note 

that within the present discrete approach the curve £l p i(q) neither touches nor crosses the border of electron-hole 
excitations, but just approaches it at the distances of the order of AE(q), and remains at this distance for q > q cr . 
To show this quantitatively let us calculate the highest excitation energy for q 3> q cr by using the expression ( pl| ) for 



the differences r2 2 (fc|| , q) — E 2 (k\\ , q) in this limit. We get from the sum rule (26) 
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In order to calculate the sum on the right-hand side we devide it in the convenient way and make the continuation, 
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After straightforward steps we get 



n(k F ,q) = E(k F ,q) 



e 2 k F 
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(30) 



noting that corrections to the energies of incoherent excitations as calculated in the previous Section [introduced in 
Eq. (|29|) through the function F(k\\,q)] do not contribute to the result ( |30| ) up to the order of w^J (q & /m?). 

To conclude, the above analysis shows that the highest branch of excitations gradually approaches the quasi con- 
tinuous incoherent electron-hole range as q approaches q cr from below. In this range it is a well-defined collective 
(plasmon) mode. For q > q cr it remains above the top of this range at the microscopic energy difference of the order 
of AE(q) or less, and as such does not have the properties of a collective mode. Simultaneously as q passes through 
q cr qualitative changes in the incoherent electron-hole range take place, as was already emphasized at the end of Sec. 
II. In Fig|| we illustrate the above discussion, which to some extent complements that from Ref. Q, with numerical 
results for a large finite system. Note that our plasmon dispersion clearly differs from those schematically presented 
in e. g. Refs. §,§. 

We close this Section with a short remark, based on the arguments given in Ref. on the calculation of Jme(q, uj) 
in the present discrete approach. In order to get a proper result for the dissipative contributions to the correlation and 
response functions, one uses the standard recipe, i. e. calculates imaginary parts of their Fourier transforms only after 
the continuation in the reciprocal q - space. This order of steps is based on the assumption that the characteristic 
energetic level spacing in the system is sufficiently small in comparison to the reciprocal time of observation of the 
system. The equivalent proposition in the classical statistical physics is that the available time for the statistical 
averaging is much shorter than the Poincare cycle time. The duration of the time of observation however becomes 
irrelevant (i. e. it can be assumed arbitrarily long), once the continuation is performed. 
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If one keeps, like in the present approach, the discrete summations throughout the calculations, the dissipative 
contributions are well-defined only after making averaging on the energy scales larger than the inherent energy level 
spacings. In particular, in our case one reproduces the correct result for Jme(q, ui) by averaging in uo on an interval not 
smaller than the energy differences AE(q) from Eq. (|^). Indeed after averaging the imaginary part of the dielectric 
function ([!]), 

2e 2 ^ 

Imefaw) = -frr g(k\\,4{8[w + E(k h qj\+5[u,-E(k h <C]}, (31) 

Q 

K|| = — k F 

over AE(q), the smallest possible energy interval consistent with the above proposition, we get 

x ,«H-AE(q)/a me 2 

Ims(q,(j) = / Jm£(q,w = — ^( ~,q). (32) 

This is the well-known result for the imaginary part of the Lindhard function in the region of the electron-hole 

quasi continuum. We note that the above averaging procedure might be particularly relevant for mesoscopic situations 
in which the cross-over from the dissipative regime to the regime without dissipation by varying the width of time (i. 
e. frequency) window becomes attainable experimentally. 



IV. CONCLUSION 



In the standard treatments with continuous wave vector the electron-hole region is reduced to a structureless 
continuum. The present analysis in which we keep systematically the wave vector discreteness leads to explicit results 
for shifts of energies of incoherent excitations with respect to the corresponding bare electron-hole excitation energies. 
Although these shifts are infinitesimal in macroscopic systems, their knowledge enables non-standard calculations of 
other physically observables like correlation energies j| and spectral functions || . The simple example of this kind is 
the calculation of the plasmon dispersion through the use of the sum rule (p6|). 

We note that the approach presented here may be straightforwardly extended to more complex macroscopic (e. g. 
multiband |^,^|) systems. Furthermore, it is obviously particularly appropriate in studies of small mesoscopic systems 
in which electron-hole excitations are characterized by essentially larger energy level spacings with respect to those in 
the macroscopic limit. Then one also encounters interesting possibility of cross-over from the dissipatively irreversible 
to reversible regime, connected with fundamental principles of thermodynamic averaging. 



APPENDIX 



In this Appendix we consider the validity of the approximation introduced by passing from Eq. (|l|) to Eq. (||). By 
this step we replace in the residual sums [represented by the function f(fc|i,q) (Eq. |l2| )] the exact value of a given 
solution f2(fc||,q) of Eq. (|l]) by the corresponding pole E(k\\,q) as defined by Eq. (|7|) . Let us start from the exact 
expression for the function F(ku, q) 



fcii — 2n/L 



5 (fcf,,q)£(fcf|,q) 
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Here the value of the zero of Eq.(]i|) is written in the form (^). The Taylor expansion of F(k\\,q) in terms of 
e(fc||,q)A£7(q) gives 



F(k {l , q) - F(k l{ , q) = Fx (fc,, , q)0(fc|| , q) A£(q) + O [[©(fc„ , q) A£(q)] 2 ] , 
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This correction can be estimated after replacing sums by integrals. We get 
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F 1 (k h q) = - 



The approximation is justified if 
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(36) 



(37) 



This condition is best fulfilled for small and large excitation energies due to the vanishing density of electron-hole 
excitations <7(fc||,q). Eq. ( |37| ) also indicates that the approximation might be less adequate in the range of wave 
numbers k in which F(k, q) = is small. The value fco for which F(ko,q) — can be determined from Eq. (p"2|). 
Assuming q <C fc F — fcp we get 



ko 1 + 17 
- — In j— 

k F l-hL 



(38) 



and ko « 0.65/c F . 

The direct insight into the validity of method follows from the comparison of analytically calculated values of 
differences ©(fcii,q) = [£l(fc||,q) — E(knq)]/ AE(q) and those obtained numerically for the mesoscopic system of 
JVw 8 - 10 6 7r/3 electrons, and L/a^ = 10, which by means of Eq. ( p2| ) gives q cr ~ 0.13fc F . It is shown in Fig.||. For 
q = O.lkp < q cr this difference raises monotonously from zero to unity, while for q = 0.3kp > q cr it shows highly non 
monotonous behavior. Finally for q = 2kp >> q cr 0(fcii,q) is close to zero in the whole range of wave numbers fcii. 
The curves in Fig. ^ clearly confirm the above estimations that the deviations of the results of analytic method from 
Sec. II from numerical calculations are negligible for small and large excitation energies. For intermediate values of 
excitation energies and for q < q cr the method is less accurate, but deviations from numerical results are still quite 
small. Finally for q >> q cr this deviation is negligible for all excitations. 
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